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O ' Abstract 

In this contribution the resuhs from various transport models on different observables 
considered as possible signals of the phase transition from hadronic matter to the 
quark-gluon plasma (QGP) - are briefly reviewed. 

^ 1 Introduction 
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The phase transition from partonic degrees of freedom (quarks and gluons) to interacting 



\0 [ hadrons is a central topic of modern high-energy physics. In order to understand the dynam- 
ics and relevant scales of this transition laboratory experiments under controlled conditions 
i/^ ■ are presently performed with ultra-relativistic nucleus-nucleus collisions. Hadronic spectra 
O ■ and relative hadron abundancies from these experiments reflect important aspects of the 
^ . dynamics in the hot and dense zone formed in the early phase of the reaction. 
Q I Estimates based on the Bjorken formula for the energy density achieved in central 

Au+Au collisions suggest that the critical energy density for the formation of a quark- 
gluon plasma (QGP) is by far exceeded during a few fm/c in the initial phase of Au+Au 
collisions at Relativistic Heavy-Ion Collider (RHIC) energies, but sufficient energy densities 
(~ 0.7-1 GeV/fm^ p]) might already be achieved at Alternating Gradient Synchrotron (AGS) 
energies of ~ 10 A-GeV [3]. More recently, lattice QCD calculations at finite temperature 
and quark chemical potential fiq |1] show a rapid increase of the thermodynamic pressure 
P with temperature above the critical temperature Tc for a phase transition to the QGP. 
The crucial question is, however, at what bombarding energies the conditions for the phase 
transition are fulfilled. Thus, it is very important to perform an 'energy scan' of different 
observables in order to find an 'anomalous' behavior that might be attributed to a phase 
transition. 

In addition to the strong interactions in the initial stage of the reaction - attributed to the 
QGP - there are also strong (pre-) hadronic interactions after/during the hadronization phase. 
Thus it becomes very important to know the impact of such (pre-) hadronic interactions on 
the final observables. The relevant information on this issue can be provided by microscopic 
transport models based on a nonequilibrium description of the nuclear dynamics [5J. 

In this contribution I present the compilation of HSD results on two of the possible 
signals of the phase transition: strangeness and charm. The HSD (Hadron-String-Dynamics) 
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Figure 1: Excitation function of the /tt^ ratio (l.h.s.) and inverse slope parameter for 
(r.h.s.) from central Au+Au (AGS and RHIC) or Pb+Pb (SPS) collisions. The solid lines 
with open squares show the results from HSD whereas the dashed lines with open triangles 
indicate the UrQMD calculations. The solid lines with stars correspond to HSD calculations 
including 'Cronin' initial state enhancement. 



transport approach [6l [7] employs hadronic and string degrees of freedom and takes into 
account the formation and multiple rescattering of hadrons; it thus dynamically describes the 
generation of pressure in the early phase - dominated by strings - and the hadronic expansion 
phase. The HSD transport approach is matched to reproduce the nucleon-nucleon, meson- 
nucleon and meson-meson cross section data in a wide kinematic range. It also provides 
a good description of particle production in p+A reactions [8] as well electroproduction of 
hadrons off nuclei |9]. In order to obtain a model independent conclusion, we also address 
the results from the UrQMD model [lOl |TT] which has similar underlying concepts as HSD 
but differs in the actual realizations. 



2 Strangeness signals of the QGP 

As has been proposed in 1982 by Rafelski and Miiller [12] the strangeness degree of freedom 
might play an important role in distinguishing hadronic and partonic dynamics. In 1999 
Gazdzicki and Gorenstein [13] - within the statistical model - have predicted experimental 
observables which should show an anomalous behaviour at the phase transition: the 'kink' - 
an enhancement of pion production in central Au+Au (Pb+Pb) collisions relative to scaled 
pp collisions; the 'horn' - a sharp maximum in the /tt^ ratio at 20 to 30 A-GeV; the 'step' 
- an approximately constant slope of K"^ spectra starting from 20 to 30 A-GeV. Indeed, such 
"anomalies" have been observed experimentally by the NA49 Collaboration [T¥| IT5]. 

In Refs. [ini [T71 [18] we have investigated the hadron production as well as transverse 
hadron spectra in nucleus-nucleus collisions from 2 A-GeV to 21.3 A-TeV within the in- 
dependent transport approaches UrQMD and HSD. The comparison to experimental data 
demonstrates that both approaches agree quite well with each other and with the experi- 
mental data on hadron production. The enhancement of pion production in central Au+Au 
(Pb+Pb) collisions relative to scaled pp collisions (the 'kink') is well described by both ap- 
proaches without involving any phase transition. However, the maximum in the K~^/tt~^ 
ratio at 20 to 30 A-GeV (the 'horn') is missed by ~ 40% [IH [IE] - cf . Fig. [1] (l.h.s.). A 
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comparison to the transverse mass spectra from pp and C+C (or Si+Si) reactions shows 
the rehabihty of the transport models for hght systems [T7j. For central Au+Au (Pb+Pb) 
collisions at bombarding energies above ~ 5 A-GeV, however, the measured m^-spectra 
have a larger inverse slope parameter than expected from the calculations. The approxi- 
mately constant slope of K"^ spectra at SPS (the 'step') is not reproduced either [T71 [18] - 
cf. Fig. [T](r.h.s.). The HSD calculations also demonstrate that the 'partonic' Cronin effect 
plays a minor role at AGS and SPS energies for the parameter T. The slope parameters 
from pp collisions (r.h.s. in Fig. [1]) are seen to increase smoothly with energy both in the 
experiment (full squares) and in the transport calculations (full lines with open circles) and 
are significantly lower than those from central Au+Au reactions for ^/s > 3.5 GeV. 

Thus the pressure generated by hadronic interactions in the transport models above ~ 
5 A- GeV is lower than observed in the experimental data. This finding suggests that the 
additional pressure - as expected from lattice QCD calculations at finite quark chemical 
potential and temperature - might be generated by strong interactions in the early pre- 
hadronic/partonic phase of central Au+Au (Pb+Pb) collisions. 

3 Charm signals of the QGP 

The microscopic HSD transport calculations (employed here) provide a suitable space-time 
geometry of the nucleus-nucleus reaction and a rather reliable estimate for the local energy 
densities achieved. The energy density e{r; t) - which is identified with the matrix element 
T°°(r; t) of the energy momentum tensor in the local rest frame at space-time (r, t) - reaches 
up to 30 GeV/fm^ in a central Au+Au coUision at ^ = 200 GeV [19] . 

According to present knowledge the charmonium production in heavy-ion collisions, i.e. 
cc pairs, occurs exclusively at the initial stage of the reaction in primary nucleon-nucleon 
collisions. The parametrizations of the total charmonium cross sections {i = Xc <//^,^') 
from NN collisions as a function of the invariant energy \/s used in this work are taken 
from [3 EDI EH [22]. We recall that (as in Refs. [211 [SI [231 |2ll 125]) the charm degrees 
of freedom in the HSD approach are treated perturbatively and that initial hard processes 
(such as cc or Drell-Yan production from A^A^ collisions) are 'precalculated' to achieve a 
scaling of the inclusive cross section with the number of projectile and target nucleons as 
Ap X At when integrating over impact parameter. For fixed impact parameter b the cc yield 
then scales with the number of binary hard collisions Nbm {cf. Fig. 8 in Ref. [22]). 

In the QGP 'threshold scenario', e.g the geometrical Glauber model of Blaizot et al. [26] 
as well as the percolation model of Satz [27], the QGP suppression '(i)' sets in rather abruptly 
as soon as the energy density exceeds a threshold value Ec, which is a free parameter. This 
version of the standard approach is motivated by the idea that the charmonium dissociation 
rate is drastically larger in a quark-gluon-plasma (QGP) than in a hadronic medium [27] . 

On the other hand, the extra suppression of charmonia in the high density phase of 
nucleus-nucleus collisions at SPS energies [281 129] has been attributed to inelastic comover 
scattering {cf. [3 [201 [21 [301 [Ml [321 [331 [311 [351 [36] and Refs. therein) assuming that the 
corresponding J/^^-hadron cross sections are in the order of a few mb [371 [Ml [39]. In these 
models 'comovers' are viewed not as asymptotic hadronic states in vacuum but rather as 
hadronic correlators (essentially of vector meson type) that might well survive at energy 
densities above 1 GeV/fm'^. Additionally, alternative absorption mechanisms might play a 
role such as gluon scattering on color dipole states as suggested in Refs. [iQl HH [121 [13] or 
charmonium dissociation in the strong color fields of overlapping strings [23] . 
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Figure 2: The J/\l/ nuclear modification factor Raa for Au + Au collisions at ^/s = 200 GeV 
as a function of the number of participants Npart in comparison to the data from for 
midrapidity (full circles) and forward rapidity (full triangles). HSD results for the QGP 
'threshold melting' scenarios are displayed in terms of the lower (green solid) lines for midra- 
pidity J/\E''s {\y\ < 0.35) and in terms of the upper (orange dashed) lines for forward rapidity 
(1.2 < y < 2.2) within different recombination scenarios (see text). The error bars on the 
theoretical results indicate the statistical uncertainty due to the finite number of events in 
the HSD calculations. Predictions for the ratio i?^^(\E'')(T^//i?^^( J/\E')(jj/^ as a function of 
the number of participants Npart are shown in the lower set of plots. The figure is taken 
from [T9] . 

The explicit treatment of initial cc production by primary nucleon-nucleon collisions and 
the implementation of the comover model - involving a single matrix element Mq fixed by 
the data at SPS energies - as well as the QGP threshold scenario in HSD are described 
in Refs. [191 EO] (see Fig. 1 of Ref. [30] for the relevant cross sections). We recall that 
the 'threshold scenario' for charmonium dissociation is implemented as follows: whenever 
the local energy density e{x) is above a threshold value Ej (where the index j stands for 
J/\E', Xc, ^')) the charmonium is fully dissociated to c + c. The default threshold energy 
densities adopted are ei = 16 GeV/fm^ for J/'if, 62 = 2 GeV/fm^ for Xc and £3 = 2 GeV/fm^ 
for Two more scenarios were implemented similarly to the 'comover suppression' and the 
'threshold melting' by adding the only additional assumption - that the comoving mesons 
(including the D-mesons) exist only at energy densities below some energy density ecut, which 
is a free parameter. We use ecut = 1 GeV/fm^, i.e. of the order of the critical energy density. 

In the following, we compare our calculations to the experimental data at the top RHIC 
energy of ^/s = 200 GeV. We recall that the experimentally measured nuclear modification 
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Au+Au, s'^^=200 GeV, comover absorption | ^ 
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Figure 3: Same as Fig. 2 for the 'comover absorption scenario' including the charmonium 
reformation channels without cut in the energy density (l.h.s.) and with a cut in the energy 
density ecut = 1 GeV/fm^ (see text for details). The figure is taken from [19]. 



factor Raa is given by, 

dN{J/^)AA/dy 
N^u-dN{J/m)Jdy' 

where dN{J /'^)aaI dy denotes the final yield of J/\E' in AA collisions, dN{J /^)ppl dy is the 
yield in elementary pp reactions and N^ou is the number of binary collisions. 

Due to the very high initial energy densities reached (corresponding to T > 2Tc), in the 
threshold melting scenario all initially created J/\E', and Xc mesons melt. However, the 
PHENIX collaboration has found that at least 20% of J/^ do survive at RHIC [44j. Thus, 
the importance of charmonium recreation is shown again. In HSD, we account for J/^ 
recreation via the DD annihilation processes as explained in detail in p[9l [30] ■ Note that in 
our approach, the cross sections of charmonium recreation in D+D —>■ J /^>+meson processes 
is fixed by detailed balance from the comover absorption cross section J /'^ +meson —>■ D+D. 
But even after both these processes are added to the threshold melting mechanism, the 
centrality dependence of the i?^^(J/\E') cannot be reproduced in the 'threshold melting' 
scenario, especially for peripheral collisions (cf. Fig. [2]). This holds for both possibilities: 
with (r.h.s. of Fig. [2]) and without (center of Fig. [2]) energy-density cut ecut, below which 
D-mesons and comovers do exist and can participate in D + D ^ J/'^ + meson reactions. 

Comover absorption scenarios give generally a correct dependence of the yield on the 
centrality. If an existence of D-mesons at energy densities above 1 GeV/fm^ is assumed, 
the amplitude of suppression of J/\Ef at mid-rapidity is also well reproduced (see the line 
for 'comover without ecut scenario in Figj3l l.h.s.). Note that this line corresponds to the 
prediction made in the HSD approach in [15]. On the other hand, the rapidity dependence 
of the comover result is wrong, both with and without ecut- If hadronic correlators exist 
only at e < ecut, comover absorption is insufficient to reproduce the J/"^ suppression even 
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at mid-rapidity (Fig. [3l r.h.s.). But its contribution to the charmonium suppression is, 
nevertheless, substantial The difference between the theoretical curves marked 'comover 
+ Ccut and the data shows the maximum possible suppression that can be attributed to a 
deconfined medium. 

We mention that there are also alternative explanations of the experimental data for 
the anomalous J/\I' suppression in A+A collisions: e.g. formation of charmonia only at the 
phase boundary as advocated by Andronic et al. [49] in the statistical hadronization model. 

4 Summary 

Summarizing this contribution, I want to point out that strange hadron production in cen- 
tral Au+Au (or Pb+Pb) collisions is quite well described in the independent transport 
approaches HSD and UrQMD [121 . The exception are the pion rapidity spectra at the high- 
est AGS energy and lower SPS energies, which are overestimated by both models. As a 
consequence the HSD and UrQMD transport approaches underestimate the experimental 
maximum ('horn') of the /tt^ ratio at ~ 20-30 A-GeV [16]. The inverse slope parameters 
T for mesons from the HSD and UrQMD transport models are practically independent of 
system size from pp up to central Pb-|-Pb collisions and show only a slight increase with colli- 
sion energy, but no 'step' in the transverse momentum slopes as suggested by Gazdzicki 
and Gorenstein [13] in 1999 and found experimentally by the NA49 Collaboration. The rapid 
increase of the inverse slope parameters of kaons for collisions of heavy nuclei (Au-|-Au) found 
experimentally in the AGS energy range, however, is not reproduced by both models (see 
Fig. [T]). Since the pion transverse mass spectra - which are hardly effected by collective 
flow - are described sufficiently well at all bombarding energies [TS], the failure has to be 
attributed to a lack of pressure. I have argued - based on lattice QCD calculations at finite 
temperature and baryon chemical potential [2113] as well as the experimental systematics 
in the chemical freeze-out parameters [48j - that this additional pressure should be generated 
in the early phase of the collision, where the 'transverse' energy densities in the transport 
approaches are higher than the critical energy densities for a phase transition (or cross-over) 
to the QGP. The interesting finding of the analysis is, that pre-hadronic degrees of freedom 
might already play a substantial role in central Au-|-Au collisions at AGS energies above 
~ 5 A-GeV. 

The formation and suppression dynamics of J/\I', Xc and ^' mesons has, furthermore, 
been studied within the HSD transport approach for Au + Au reactions from FAIR to 
top RHIC energies of ^/s = 200 GeV. It is found that both the 'comover absorption' and 
'threshold melting' concepts fail severely at RHIC energies [T9] whereas both models perform 
quite well at SPS energies. The failure of the 'hadronic comover absorption' model goes in 
line with its underestimation of the collective fiow V2 of leptons from open charm decay 
|46j . This suggests that the dynamics of c, c quarks at RHIC energies are dominated by 
strong pre-hadronic/partonic interactions of charmonia with the medium in a strong QGP 
(sQGP), which cannot be modeled by 'hadronic' scattering or described appropriately by 
color screening alone. 

The evidence for creating a 'new state of matter', the sQGP, in Au+Au collisions at RHIC 
is overwhelming and is additionally supported by the strong suppression of high pt hadrons 
and 'far-side' jets - both observations being insufficiently described by string/hadronic models 
[501I5T] - as well as the quark-number scaling of elliptic fiow V2{pt)- The question now reads: 
what are the properties of the sQGP? 
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